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We examine the role of order-parameter suppression in the development of low-energy peaks 
(i.e., resonances) in the tunneling density of states near a non-magnetic impurity in a d-wave su- 
perconductor. Without order-parameter suppression, the zero-energy resonance appears only in 
the unitary (i.e., strong impurity) limit. However, suppression makes the resonance appear even 
when the impurity is much weaker. To model this situation, we make the physical hypothesis that 
the order parameter is reduced whenever one electron of a Cooper pair encounters the impurity, a 
hypothesis that retains the exact solvability of the problem. In this way, we determine that suppres- 
sion of the order parameter drives the effective strength of the impurity towards the unitary limit. 
We determine the order-parameter reduction variationally, and show that the ratios between the 
main energy scales — the band width and superconducting gap — strongly affect this reduction and, 
in consequence, the position and width of the resonance. 



I. INTRODUCTION 

The role played by non-magnetic impurities in high-temperature superconductors (HTSCs) represents an important 
element of the subject of high-temperature superconductivity. In contrast with the case of conventional superconduc- 
tors, in which the s-wave symmetry of the order parameter tends to weaken the effect of such impurities (c.f. Anderson's 
theoremEl), the HTSC materials display rich and interesting sensitivity to the amount of disorder, even at low disorder- 
concentrations. Indeed, many physical properties are affected at low temperatures and frequencies, the most direct 
example being the appearance of a nonzero density of states (DOS) at the Fermi levelo. One of the reasons for 
this sensitivity lies in the properties of individual non-magnetic impurities in a d-wave superconducting host, such 
impurities giving rise to resonant quasiparticle states at-.sub-gap energies. The occurrence of these states was pre- 
dicted theoretically by Balatsky, Salkola and co-workersau, and is consistent with recent experimental observations 
by Yazdani et al.tj These states are localized near impurity sites, and have finite life-times, due to the existence of 
bulk quasiparticle states into which they may decay. Upon neglecting the changes in the order parameter induced 
by an impurity, it was found that, as the strength of the impurity increases, the resonances move towards the Fermi 
level and their widths decrease; only in the unitary (i.e., infinitely-strong impurity) limit do the resonances reach the 
Fermi level and become infinitely sharp. The role of the order-parameter changes (i.e., suppression) has been analyzed 
by several groups (see, e.g., Refs. [ In the present Paper we focus on one particular effect of order-parameter 

suppression which has not been reported previously: we show that the suppression of the order parameter drives (i.e., 
renormalizes) the effective strength of the impurity towards the unitary limit. Further, we argue that this renormal- 
ization may be appreciable in the HTSC materials, and may be important for the development of a more complete 
understanding of the low-temperature behavior of the cuprates. Indeed, as was argued in Ref. [ px| , the observed 
low-temperature behavior of cuprates is inconsistent with the relatively small suppression of the critical temperature, 
' unless the impurities are (or at least behave as if they are) in the unitary limit. 

II. MODEL OF ORDER-PARAMETER SUPPRESSION NEAR AN IMPURITY 

_ _ i 

Consider a point-like impurity in a two-dimensional d-wave superconductor. As a fully self-consistent treatment of 
the order-parameter suppression is out of reach, we shall proceedrby exploring a physically motivated hypothesis for 
the functional form of this suppression. It is commonly assumeduH that this suppression takes the form 

5A( r y)ocf( T -^-)A (r-r'), (1) 
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where / gives the spatial shape of the suppression, Ao(r — r') is the bulk value of the d-wave order parameter, and the 
impurity is located at the origin. This form includes only the d-wave pairing channel, and therefore is very convenient. 
We show, however, that Eq. (jl|) is meaningful only for smooth / (varying on length scales much longer than the Fermi 
wave length kp 1 ), and does not describe the physical situation at short distances from the impurity. This short- 
distance behavior is important, as it affects the formation of the scattering resonances. As an extreme example, 
consider the setting of tight-binding electrons moving on a two-dimensional square lattice with on-site repulsion and 
nearest- neighbor attraction (i.e., the simplest situation for d-wave superconductivity). If we locate an impurity at the 
origin and wish to suppress the order parameter in the four bonds connecting the origin to its nearest neighbors, we 
arrive at the following functional form of the suppression: 

SA(r, r') = a( 6{r) + 6{r') ) A (r - /) , (2) 

where a is the amplitude of the suppression which has the dimensionality of a volume. This form may also be used 
for other cases, inasmuch as it encodes the idea that the order parameter is altered whenever one of the electrons in 
the Cooper pair encounters the impurity. The Fourier transform of Eq. (|^) reads 

SA(k,k') = J drdr' 6A(r,r') exp{ikr + ik'r') = a (A (k) + A (-fc')) • (3) 

As usual for superconductivity, the most important regime is the one in which both k and k' are close to the Fermi 
surface. Thus, the assumption (|^) may be relaxed in favor of (Q) near the Fermi surface. The form (||) includes pairing 
channels other than d-wavc. Let us, e.g., adopt the tight-binding shape of the order parameter: Ao(fc) = A0d(fc), 
where 



0d/ s (fc) = cos(ap x ) T cos(ap y ) , (4) 

in which a is the lattice constant, and the subscripts d/s stand for the d and extended-s channels. By introducing the 
total [q = k + k'] and the relative [p = (k — k')/2] momenta of a Cooper pair we arrive at 

5\, q = a A (&i(p)&(g/2) + Mp)M<i/2)) ■ (5) 

Recent numerical self-consistent simulations! yield an order-parameter suppression having a form similar to (||) , with 
the d-wave contribution having the form-factor of the s-wave symmetry and vice versa. Note that at small values 
of q only the d-wave contribution survives, which reconciles Eqs. (||) and (0). Further support for the choice of 
the functional form (|J) may be provided by examining the imbalance in the self-consistent equation for the order 
parameter with the impurity present and the order parameter unchanged (i.e., the model considered in Ref. ^). The 
imbalance means that the gradient of the free energy in the function space of A(fc, k') is nonzero and the "direction" 
of this gradient gives the functional form of the linear response of the order parameter SA to the presence of the 
impurity. We find that the "direction" of the imbalance is close to (p|) , and also has the same symmetry as M) . 



III. DENSITY OF STATES NEAR THE IMPURITY 



First, we investigate the local DOS near the impurity. We employ the standard T-matrix technique. The system is 
governed by the Hamiltonian 

H{r, r') = H (r, r') + E(r, r') = H (r, r') + U5(r)S(r')a z - 6A(r, r')a x , (6) 

where Hq is the Bogolubov-de Gennes kernel, which describes the unperturbed d-wave superconductor, and the term 
with the coefficient U represents the potential scatterer (i.e., impurity), hats denote 2x2 matrices in the Nambu 
space. 

The Dyson equation for the full Matsubara Green function G(r, r', iu> n ) = (ioj n — H(r,r'j) reads 



G(r,r')=G (r,r , ) + J dx ± dx 2 G (r,x 1 ) E(x u x 2 ) G(x 2 ,r') (7a) 
= G (r,r')+ dx 1 dx 2 Go(r,x 1 )f(x 1 ,X2)G (x2,r'). (7b) 
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Our aim is to find the T-matrix T(xi,x 2 ), which would allow us to calculate, via Eq. (7b), the full Green function 
and, therefore, the DOS. The Dyson equation for the T-matrix reads: 



f(xi,x 2 ) = t(x 1 ,x 2 ) + / dy 1 dy 2 t,(xi,yi)Go(yi,y2)f(y 2 ,X2) 



(8) 



This model is exactly solvable, owing to the fact that the Fourier transform of the self-energy £(fc, fc') is a degenerate 
kernel (i.e., a sum of factorized functions of k and k'): 

(9) 



E(k, k') = Ua z ~a (A (k) + A (-k'))a x . 
By using this degeneracy property we may rewrite Eq. (||) as 

T(fc, fc') = Ua z - a(Ao(fc) + A (fc')) &x + U&, F{k') -aa x H(k') - a A (k) a x F(k') , 

F{k')=V- x Y,Go{k")f{k",k') , 
fc" 

H{k') = V- l Y,Mk")G Q {k")f{k",k') , 



(10a) 
(10b) 

(10c) 



the remaining task being to determine the as-yet unknown ma trix- valued functions F and H. Here and subsequently, 
we use the symmetry Ao(— fc) = Ao(fc). By multiplying Eq. (10a) first by Go(fc) and then by Ao(fc)Go(fc) from the 
left, and integrating over fc, we obtain a system of linear equations for F and H: 

(UPa z -aLa x -l) F(k) ~aPa x H(k) = -UPa z + a L a x + a P A (fc) a x , (11a) 
(UL& Z -aM a x ) F(k) - (ata x + l) H{k) = -ULa z + a M & x + a L A (k) a x , (lib) 

where P = V' 1 J2 k G (fc), L = V' 1 J2 k A (k) G (fc), and M = V' 1 J2 k A 2 (k) G (fc). 

Thus far, we have not made use of the d-wave symmetry of the order parameter, and we have not made any 
approximation beyond the mean-field (Bogolubov-de Gennes) approximation. To determine the matrices P, L and M 
we use the d-wave character of the order parameter to eliminate integrals of the odd powers of Ao(fc), and we assume 
the presence of particle-hole symmetry to eliminate integrals of odd powers of the single-electron energy e(fc). Taking 
for the unperturbed Green function 



Go(fc) 



A%(k) + e*(k)+u,Z 



-iu> n -e(k) -A (fc) 
-Ao(fc) -iuj n + e(k) 



we obtain P(iuj n ) — P(icu n )(— 1), L(iuj n ) — L(iuj n )(—a x ), and Al(iuj n ) — iLu n L(iuj n )(—l), where 

P(iw n ) = V- 1 



A2(fc)+ e 2(fc)+C,2 



L(iLo n ) = V^ 1 ^2 



Ag(fc) 



-A§(fc) + e 2( fc)+ ^ 



(12) 

(13a) 
(13b) 



The system ( 11a , lib ) consists of a pair of 4-dimensional systems of linear equations having common coefficient 
matrices and distinct inhomogeneous terms. These two systems may be analytically solved,and thus the exact T- 
matrix may be rebuilt using Eq. ( 10a ). Then the exact Green function can be built by using (|7b|), and hence the local 
DOS may be calculated via 



(r,E) = --ImG 1A (r,r,E + iS)\ S - 



(14) 



Before we present exact results for the DOS, it may be noted that the singular behavior (i.e., resonanc es) of th e 
T-matrix at sub-gap energies (for E < A) origi nate only from zeros of the determinant D(E) of the system (11a, lib). 
This is so because the right hand side of Eqs. (11a, lib) has at most branch-cuts at the sub-gap energies. Thus, it is 
instructive to write down this determinant: 



D(E) = D 1 {E)D 2 (E) , 
D lt2 {E) = 1 - 2aL(E) + a 2 L 2 (E) - a 2 EL(E)P{E) ± UP(E) 



(15) 
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If a — (i.e., no order-parameter suppression) the characteristic equation D(E) = is identical to the equation for the 
poles of the T-matrix obtained in Ref. [ Thus, we may conclude that the role of the order-parameter suppression 
in the present model is to modify (i.e., shift) the resonance found in Ref. [ 0], rather than to add a new resonance (as 
results from the model proposed in Ref. [ m) . 

We now calculate the DOS precisely at the impurity. We obtain 



v(r, E)\r=o = v Q {E) Im 

TT 



P(E)(D 1 (E)-l 
Di(E) 



= — Im 


\ P(E) ) 


TT 


[Di{E)\ 



(16) 



where Vq(E) = (tt r ) Im P(E) is the unperturbed DOS in the d-wave superconductor. For a — 0, the behavior of 
v(r,E) for r = was described in Ref. [ |J. There, it was shown that a resonant peak appears at negative energies 
when U = Uvp becomes comparable to 1 (up being the DOS at the Fermi surface in the normal state). The peak 
moves toward zero energy and becomes narrower and weaker as U grows. In the unitary limit (i.e., for U — > oo) the 
peak disappears. This disappearance simply means that an infinitely strong impurity repels all the electronic density 
from itself. The resonance is still there, and to analyze it, one should calculate the tunneling density of states in the 
vicinity of the impurity (see Ref. [ Q). There, four maxima appear along the lobes of the d-wave order parameter at a 
distance of the order of the Fermi wave length Xp from the impurity site. Moreover, a second resonance, corresponding 
to a singularity of the subdetcrminant D%, shows up in the vicinity of the impurity. The width of this second resonance 
is exactly equal to the width of the first, and the positions of the two are symmetric with respect to the Fermi energy. 
However, the spatial density distributions of the two resonances differ from one another. 

The new effect that we report here is that if one fixes the impurity potential U and allows the order-parameter 
suppression a to grow instead, the DOS behaves similarly to the scenario outlined above. Specifically, the resonances 
move toward zero energy, become sharper, and a maximum in the DOS is fou nd al ong the lobes of Ao(fc) at a distance 
of the order of from the impurity. Indeed, for E — > 0, we obtain from Eq. ( |l3tf ) that L(E) — > 2Avp [we have taken 
Ao(fc) = Acos(20fc)]. In this regime we can approximate D\^(E) as (1 — 2a)' 2 ± UP(E), where a = aAvp. We find 
that the positions (±Qo) and the width (r) of the resonances may be now determined using the formulas of Ref. [ Q: 



n 



A 



2U eS ]n(8U, 



r = 



ttQo 



cir ) 



21n(8C/ cff ) 



(17) 



in which the original strength of the impurity is substituted by a renormalized one, viz., U c fr = U/(l — 2a) 2 . Now, 
unitary behavior is achieved if U e s S> 1, and a strong renormalization of the bare strength of the impurity occurs 
if it is possible for a to be close to 1/2. In the next section we will argue that such a regime may be realistic in 
the HTSC materials. In the meantime, let us assume that strong renormalization has occurred, i.e., that the bare 
strength of the impurity U was not large enough to cause the unitary behavior but that U e s is. Then, we may ask 
the question: Is there any difference in the spatial distribution of the resonant-state density between this case and the 
case when [/>1? We find that the only difference is precisely at the impurity site: there, the electrons "know" that 
the impurity is not so strong, and therefore the DOS is less suppressed than in the "true" unitary limit (i.e., U ^> 1). 
Farther from the impurity, however, the two cases are indistinguishable. The scaling relation between these two cases 
may be expressed as 



i/(r, E, U , a) 



v{r, E,U cS ,a) | a = 
(1-25)" 



for r = 0, E « 0; 



(18) 



f(r, E, U cS , a) | s=0 , for \r\ > X F , E « 0. 



It is not clear that such a behavior could be detectable experimentally because even when STM tips are precisely 
above impurities tunneling occurs over some neighborhood of the impurity. However, this is at least consistent with 
the experimental results reported in Ref. [ |j| . 



IV. AMPLITUDE OF THE ORDER-PARAMETER SUPPRESSION 



Our next step is to estimate the amplitude of the order-parameter suppression a. In principle, one may envision 
two different scenarios. In the first, the electron-electron interaction is unchanged by the presence of the impurity, 
and, in the d-wave case, the suppression of the order-parameter is only due to the pair-breaking effect of the impurity. 
In the second scenario, the electron-electron interaction is itself suppressed near the impurity, thus furthering the 
suppression of A(r, r'). Let us make a very crude estimate for the second scenario. We again exploit the tight-binding 
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model, and assume that the order parameter is zero in the four bonds connecting the impurity to its neighbors but 
unchanged elsewhere. This regime would be reasonable for purely electronic mechanisms of superconductivity, as the 
local electronic structure is completely altered by the impurity. Then a ~ a d and a — a Av m A/2t, where 2t is the 
band- width and we have used as an estimate for the density of states v « l/2ta d . If we take into account the fact 
that in the definition of a the DOS at the Fermi level Vp should be used, whereas the lattice constant a is naturally 
connected to the DOS averaged over the whole band (i.e., D), we arrive at a more refined estimate: a s=s (A/2t)(i'p/D). 
We see that in conventional superconductors a is always small, and thus essentially no renormalization can happen. 
In the HTSC materials, however, A/2t can be of the order 0.1, and the proximity of the van Hove singularity makes 
the factor vpjv important. Thus, in this case a strong renormalization situation cannot be ruled out. 

Next, we show that, even if the electron-electron interaction is unchanged near the impurity, the pair-breaking 
process creates a suppression of the order parameter to the value estimated above [a « {A/2t){vp/v) ] when U ~ 1. 
We establish this variationally, i.e., we minimize the free energy of the system with respect to a. To calculate this 
free energy one has to know the form of the electron-electron interaction responsible for the superconductivity of the 
system. 

The most general form of the pairing interaction may be written as 



VN/2,T °-p+q/2,! L -p'+ g '/2,l V+9'/2,T 



(19) 



P,P' ,1,1' 



Note that we have introduced the Kronecker symbol Sq, q >, corresponding to the conservation of the total momentum, 
in order to emphasize the matrix structure of this interaction. This matrix has a two-fold index (p, q) corresponding 
to the relative and the total momenta of the electrons in a Cooper pair, respectively. To avoid confusion, we shall 
use the letters p and q for the relative and the total momenta of the electrons in a Cooper pair, -reserving the letter k 
for momenta of individual electrons. The standard Hubbard-Stratonovich decoupling procedureEU yields the following 
effective action for the order parameter: 



S = Si + S2 — / dr 



V 3 



p,p',i,q' 



-TrlnG" 1 



(20) 



where G 1 (r, t; r', t') = iS(t—t')S(r—r') id/dt' — 5{t — t') H(r, r'), and g(p,p'; q) stands for the inverse of the g(p,p'; q) 
matrix in p space [the inversion in q space is trivially performed in Eq. (EQ)]: 



(21) 



The electronic free energy F in the mean-field approximation is thus given by 

S 



F[A p , q ] = ^ = F 1+ F 2 = -±Y, 9(p,p';q) \, q ^. q ~ ^TrlnG" 1 



V 3 



p,p',g 

We now minimize (|2^) with respect to a. The second term of ( p2|) is treated readily; indeed 

1 



da 2 da \ /3 ^ 



-TtG—G^ 1 

(3 da 



(3V 2 



k,k' 



tr G(fc,fc',z Wn )(A (fc) + A (fc'))^ 



(22) 



(23) 



where the symbol tr stands for the trace in the Nambu space only. The Green function G is known exactly: G(k, k') 
V5 k ,-k' G (k) + G (k) f(k, k') G {k'). Therefore, we can rewrite Eq. (|| as 



(3V 



k,u) n 



13V 2 



Y I A o(fc) + A (fc')] tr G (k,iuj n )f(k,k',iu; n )GQ(k',iuj n )a x 



k,k' 



As for the first term of Eq. (p2|), we recall that 

A p , g = V6 qfl A (p) - a (Ao(p + q/2) + A {p - q/2)) 



(24) 
(25) 
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and therefore F\ is a quadratic polynomial in a. The derivative of F\ with respect to a contains a term independent of 
a and a term linear in a. It is straightforward to verify that the term independent of a in dF\ /da, when combined with 
the first term in Eq. (pi|), cancel, as together they constitute the BCS self-consistency equation for the unperturbed 
order-parameter. (The appearance of this equation was to be expected as it emerges from the condition that the free 
energy be minimal when no impurity is present.) We denote the parts of F\ and Fi remaining after the cancellation 



as F\ and F^' , respectively, and the variational condition for a now reads: 



.(2 



dF 
da 



(2) 



± F 

da 1 
1 



(2) 



da ' 



(2) 



13V 2 ^ 

k,k .uj 7 



2A (fc)tr Go(k, iuj n ) T(k, k', iu n ) Go(k', iw n ) cr 2 



r(2) 



1 



p,p',q 



^ g(p,p';q) 5A Ptq SA* p , iq 



To estimate F, 



(2) 



we choose the standard form of the interaction, which has no g-dependence: 

g(p,p'\ q) = -9d <t>d(p) <t>d{p) - g s Mp) Mp) ■ 

We include here both d and s channel interactions. The part of F\ quadratic in a reads 

, ( 2) i 1 1 MpHMp') . i *,(p)*,(pT 



9d 



(26a) 
(26b) 

(26c) 

(27) 
(28) 



where N, 



To proceed further we need to assume some particular form of the wave functions 4>d/s- The simplest choice is the 
tight-binding one [Eqs. (Q) and (EM, for which F^ is readily calculated: 



F (2) = « 



(9^+ 9s 1 ) 



(29) 



where we have used the identity V 1 J2 q <p 2 l / s (q/'^) — l/a d - The appearance of the volume of the lattice cell a d in 
Eq. ( p9| ) introduces the band-width energy-scale 2t, via v m l/2ta d . Defining, as usual, the pair of dimensionless 
coupling constants <?(j/ s = vp gd/s; an d differentiating Eq. ( p9| ) with respect to a, we obtain 



d_ 

da 



F (2) = 



2aA 2 



{9a 1 



2«(ffd" 1 +57 1 ) 



2W 



A 2 u F 



(30) 



Finally, to find a we must evaluate dF^ /da, as given by Eq. (26b). Although all the components of ( |26b| ) are 
known a naly tically and the integrals over k and kl can be expressed via the functions P{iuj n ) and L(iu> n ) defined 



by Eqs. (13a) and (13b), the remaining sum over the Matsubara frequencies must be carried out numerically. An 
analytical result is obtained only in the "true" unitary limit U — * oo: 

d_ 

da 



(31) 



where the last equation is obtained from the self-consistency condition for the unperturbed order parameter (without 
impurities). ^From Eqs. (^0|) and (|3l| ) we obtain 



25d' 



A v F 



gs 1 2t v 



A v F 
2t ~V 



(32) 



This result would still be meaningless if the asymptotic value (|32|) were achieved only for U oo, as the renormal- 
ization has no effect in the "true" unitary limit. To check how fast this asymptotic value is achieved, we have solved 



Eq. (26a) numerically for different choices of the system parameters. We observe that the asymptotic value is always 
reached already for U f» 1 (see, e.g., Fig. |l|). Thus, impurities having "mild" values of U may be renormalized to the 
unitary limit. We were unable, however, to approach numerically the regime awl/2 without employing the van Hove 
singularity (i.e., for vp/v m 1). This is because the ratio A/2t would have to become of order 1, which is inconsistent 
with the BCS approximation. On the other hand, exploring the van Hove scenario would demand calculations with 
realistic band structures (see, e.g., Ref. [ |J), the task we leave for the future. 
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V. CONCLUSIONS 



To conclude, we have investigated the role of the order-parameter suppression in the development of resonant 
scattering states around impurities in d-wave superconductors. We show that the suppression of the order parameter 
drives the effective strength of the impurity towards the unitary limit. This effect may be relevant in the HTSC 
materials, due to the relatively large value of the ratio between the superconducting gap A and the band-width and 
due to the possibility of a van Hove singularity in the DOS. The electronic DOS around a renormalized impurity is 
indistinguishable from the DOS around a "truly" unitary impurity, except precisely at the impurity site. Further 
calculations involving the strong-coupling regime and the effects of the real band-structures are needed to establish 
the feasibility of the strong renormalization of the strength of impurities. 
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FIG. 1. Results of the numerical simulations for a(U). A/2t = 0.1, vf /v = 1, the coupling constant §d « 0.8 is determined 
self-consistently, and it is assumed that d s = dd- 
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